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ABSTRACT 

Measurements have been made of the propulsive effect of supersonic combustion ramjets 
incorporated into a simple axisymmetric model in a free piston shock tunnel. The 
nominal Mach number was 6, and the stagnation enthalpy varied from 2.8 MJ kg' 1 to 
8.5 MJ kg’ 1 . A mixture of 13% silane and 87% hydrogen was used as fuel, and 
experiments were conducted at equivalence ratios up to approximately 0.8. The 
measurements involved the axial force on the model, and were made using a stress wave 
force balance, which is a recently developed technique for measuring forces in shock 
tunnels. A net thrust was experienced up to a stagnation enthalpy of 3.7 MJ kg' 1 , but as 
the stagnation enthalpy increased, an increasing net drag was recorded. Pitot and static 
pressure measurements showed that the combustion was supersonic. 


The results were found to compare satisfactorily with predictions based on established 
theoretical models, used with some simplifying approximations. The rapid reduction of 
net thrust with increasing stagnation enthalpy was seen to arise from increasing 
precombustion temperature, showing the need to control this variable if thrust 





1 


\ 

performance was to be maintained over a range of stagnation enthalpies. Both the 

1 1 

| inviscid and viscous drag were seen to be relatively insensitive to stagnation enthalpy, 

ii 

■ with the combustion chambers making a particularly significant contribution to drag. The 
! maximum fuel specific impulse achieved in the experiments was only 175 sec., but the 

theory indicates that there is considerable scope for improvement on this through 

/. 

aerodynamic design. /// 

_ (/' 

1. INTRODUCTION 

The supersonic combustion ramjet, or "scramjet", has existed as a concept for some 
decades (e.g. Fern 1964, Swithenbank 1967, Jones and Huber 1978, Billig 1993). 
Essentially, it generates a propulsive effect through heat addition, by combustion, at 
supersonic speeds. In principle, it is able to operate at any flight speed, but a practical 
limit seems to exist at 4 or 5 km s' 1 . To explore the potential of the scramjet, the ability 
to conduct routine experimentation at such speeds is essential. 

For speeds in excess of 2 km s\ power and materials limitations have forced ground 
facilities to adopt the principle of impulse operation, with test times of only a few 
milliseconds. The most successful facility of this type has proven to be the shock tunnel, 
and indeed, substantial scramjet experimentation has been done in shock tunnels (e.g. 
Stalker and Morgan 1984). 

However, the scramjet is a propulsive device, and the ultimate test of a propulsive device 
is its ability to generate thrust when installed in a flight vehicle. Therefore 
experimentation in which thrust or drag is measured on vehicle models in a shock tunnel 
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can be expected to be an important part of scramjet research.. 


Methods of measuring forces on a shock tunnel model have been reported by previous 

* 

investigators (Bernstein 1982, lessen and Gronig, 1989, Naumann et al 1993), but these 
are only effective with test models which are relatively small with respect to the tunnel 
size. For a scramjet model to produce thrust, vigorous combustion in the combustion 
chambers is required and, with the combinations of pressure and Mach number which can 
be achieved in a shock tunnel, a not insubstantial length is required for the combustion 
chambers. Since the combustion chambers must be integrated into a complete model it is, 
therefore, very difficult to retain a model size which is small in relation to the tunnel 
size. 

This paper reports measurements of the effect of stagnation enthalpy on the thrust and 
drag of a quasi-axisymmetric scramjet model. It is of particular relevance to the possible 
use of scramjets to propel acceleration vehicles, where one of the effects of the 
continually changing flight velocity will be a changing stagnation enthalpy. The model 
incorporated intakes, combustion chambers and thrust nozzles, and was therefore 
aerodynamically complete. Having regard to the difficulty noted in the previous 
paragraph, the measurements were accomplished by employing a stress wave force 
balance (Sanderson and Simmons 1991), which will be explained in more detail below, 
and which does not require a model which is small in relation to the tunnel size. A 
positive resultant thrust was measured - i.e. the total thrust generated was greater than the 
drag of the complete configuration - though only over a limited range of test section 
velocities. 


3 


The measurements are also compared with theoretical predictions. Though intakes, 
combustion chambers and thrust nozzles have all been the subject of research separately, 
an experimental study involving all three elements combined into an integrated 

o 

configuration has not been reported, at least, in the open literature. Therefore the 
comparison between theory and experiment is of considerable interest. The model was 
designed to facilitate this comparison, and so an outstanding propulsive performance was 
neither expected nor achieved. 

The shock tunnel in which the experiments were conducted is briefly described first, 
paying attention to limits on the test time and the test section flow produced. Then the 
main features of the model are described, followed by a short discussion of the principle 
of the stress wave force balance. The application of the balance to this problem is 
outlined, and some typical data recordings are presented. Following this, theoretical 
predictions of the effect of stagnation enthalpy on the thrust are developed. These 
predictions are then compared with measurements of the thrust, the comparison being 
assisted by some pressure measurements in one of the combustion ducts of the model. 

2. THE SHOCK TUNNEL 

The experiments were done in the free piston shock tunnel T4 at The University of 
Queensland, as shown in fig. 1. This consisted of a compression tube, 27 in long and 
228 mm in diameter, in which the temperature and pressure of the driver gas was raised 
by a free piston compression process immediately prior to diaphragm rupture. The shock 
tube was 10 m long and 76 mm in diameter, and for these tests, supplied test gas to a 
contoured hypersonic nozzle which expanded from a throat diameter of 25 mm to a 
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diameter of 263 mm at the test section. 


The nozzle stagnation pressure was monitored during each test by a PCB quartz 
piezoelectric pressure transducer located 60 mm upstream of the downstream end of the 
shock tube, and a typical record is shown in fig. 2(a). It can be seen that after the initial 
peak at shock reflection, the pressure remained constant within 20% for approximately 2 
millisec. A mixture of Argon with Helium was used as the driver gas, and the ratio of 
the two was varied as the stagnation enthalpy was varied to ensure that a similar period of 
approximately constant nozzle stagnation pressure was maintained throughout the tests. 

Measurements with a time-of-flight mass spectrometer on the nozzle axis were used to 
determine when the test section become contaminated by driver gas, and results are 
shown in fig. 2(b)*. It can be seen that, although the time to contamination exhibits 
substantial shot to shot variation at a given stagnation enthalpy, a period of 1 millisecond 
or more of flow with a contamination level of less than 10% exists up to stagnation 
enthalpies of 10 MJ kg' 1 . This is in excess of the stagnation enthalpies used in the 
experiments. 

The test section conditions were calculated using a non-equilibrium nozzle code (Lordi et 
al 1966) for a steady expansion in a hypersonic nozzle from known reservoir conditions 
to a given Pitot pressure level. The reservoir conditions were determined by first 
obtaining the pressure and stagnation enthalpy after shock reflection from measurements 
of the shock speed and shock tube filling pressure. It is then assumed that an isentropic 

*The authors would like to thank K. Skinner for the use of this data. 
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expansion to the measured nozzle reservoir pressure occurs. This pressure varied from 
34 MPa to 40 MPa, with an underlying trend to increase with increasing stagnation 
enthalpy superimposed on a random shot to shot variation. This is reflected in the Pitot 
pressure measurements shown in fig. 3(a). The ratio of Pitot pressure to nozzle reservoir 
pressure was (14.1 ± 0.6) x 10" 3 , and showed no variation with stagnation enthalpy so 
the scatter in Pitot pressure is due largely to variation in nozzle reservoir pressure. Using 
the value of 14.1 x 10' 3 for the nozzle calculations, the test section conditions shown in 
fig. 3(b) were obtained. 

The test section flow was calibrated for uniformity and flow divergence by traversing 
with a Pitot rake. A typical result is shown in fig. 4.’*. Two traverses, taken at 
100 mm and 300 mm from the nozzle exit respectively, are shown in the upper half of 
the figure. The position of the model is shown in the lower half. The cone of uniform 
flow starts well upstream of the model, and the traverses show that a reasonably uniform 
test core, without observable flow divergence, persists to a distance well downstream of 
the model. The location of a Pitot probe, which was used to monitor the Pitot pressure 
during tests on the model, is also shown in the bottom half of the figure and yielded the 
results of fig 3(a). 

3. THE MODEL 

The scramjet model is shown in fig. 5 with half of the cowl removed. It consists of an 
axisymmetric centrebody, with six combustion chambers and associated intakes 
symmetrically arranged about its periphery, and separated from each other by six 

“The authors would like to thank S. Overton for the use of this data. 
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splitters. The centrebody has a conical forebody with an 9° half angle, and a conical 
afterbody of 10° half angle. The combustion chamber intakes consist of compression 
ramps formed by the splitters which deflect the flow in the shock layer on the conical 
forebody through an angle of 8° parallel to its surface, with the leading edge of the cowl 
swept at an angle chosen to trap the nearly two-dimensional shocks so formed, and 
prevent them spilling into the freestream. The internal surface of the cowl is parallel to 
the forebody surface, and the external chamfer on the cowl leading edges makes an 
included angle of 8.5° in the freestream direction. The leading edge of the compression 
ramps is 10 mm wide, measured in the radial direction. 

After the 'flow in the forebody shock layer has passed through the shock formed by one of 
the compression ramps, and its reflection, it is then turned to become parallel to the 
freestream, and enters one of the combustion chambers. These are 59 mm long, and each 
has a cross-section which is an arc of an annulus of major diameter 67 mm and minor 
diameter 51 mm, subtending an angle of 23° at the centreline. The fuel injection orifices 
were located at the beginning of the cylindrical part of the centrebody, and were therefore 
at the upstream end of the combustion chambers. There was one orifice to each 
combustion chamber. They were 2 mm in diameter, angled at 30° to the centreline, and 
injected in the downstream direction. 

At the downstream end of the combustion chambers, the sidewalls deflect through an 
angle of 12°, and the conical afterbody begins, to make an array of thrust generating 
surfaces. These are terminated when the 12° surfaces meet, and when the afterbody 
diameter becomes 25 mm, at which point the centrebody becomes cylindrical again. The 
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trailing edge of the cowl also makes an internal angle of 7° with the freestream direction, 
and this is expected to add t small increment of thrust. 

The model as it is mounted in the test section is shown in fig. 6 . The model is attached 
via the valve system to the fuel tank, which is itself attached to the stress wave bar of the 
stress wave balance. The stress wave bar is suspended from ports in the tunnel dump 
tank by two fine wires. The valve assembly, the fuel tank and the stress wave bar are 
surrounded by an aerodynamic shroud, and the only contact of the model assembly with 
the shroud is through a flexible "O" ring seal on the parallel section downstream of the 
conical afterbody of the model. A streamwise section of the model assembly is shown in 
fig. 6 . It can be seen that a 10 mm hole is drilled along the model axis to carry the 
fuel, which is distributed to the injectors as shown. 

A schematic depicting the operation of the valve is shown in fig. 6 . The solenoid valve, 
S, is normally closed, and the fuel tank is filled before a test. The cavity G is maintained 
at the fuel pressure, so the fuel valve F is closed. When the tunnel is fired, its recoil 
closes a switch which activates the solenoid valve, causing it to open, and allow the fuel 
to flow out of cavity G to vacuum. The pressure in cavity G falls, the fuel valve F 
opens, and fuel flows through it to the fuel pipe. As it does so, its pressure is monitored 
by a PCB type 11 1A quartz piezoelectric pressure transducer. When the pulse of test gas 
has passed, the solenoid valve closes again, the fuel continues to flow into cavity G 
through the bleed orifice O, the pressure in the cavity builds up, and valve F closes. 

The centrebody and splitters of the model are made of aluminium alloy, while the cowl 
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was made of stainless steel. The fuel tank and valve assembly is made of stainless steel. 


4. THE STRESS WAVE FORCE BALANCE 

It will be observed that the length of the model is not small with respect to the tunnel 
size. Also it is noted by Sanderson & Simmons (1991) that stress waves, induced in the 
model by the flow starting process, propagate along the model with speeds of the order of 
3000 to 5000 m s' 1 , which is of the same magnitude as the flow speed. Therefore, the 
many stress wave reflections which must occur before the model can be said to behave as 
a rigid one will correspond to a large number of flow passes over the model, requiring 
a time which normally exceeds the test time. 

Fortunately, the stress waves which are the source of this difficulty can themselves be 
exploited to measure the force on the model. If strain gauges are placed on the model 
support sting adjacent to the base of the model, and the time history of stress waves 
passing into the sting is recorded, this recording can be deconvoluted to yield the force on 
the model. The simplest case of deconvolution occurs when the model is simply an 
extension of the stress bar and the force is applied as a step change in time at its upstream 
end, in which case the strain gauges record the passage of a stress wave which directly 
represents the applied load. This, in fact, is identical to the bar gauge which is used for 
measuring transient pressures. However as the models become more complicated, so 
does the deconvolution. 

When Sanderson and Simmons first proposed this technique, they applied it to 
measurement of drag on a short cone. In that case the internal stress wave effects were 
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unimportant, but subsequent measurements on a 425 mm long cone, (Tuttle, 1990) where 
they were important, showed that they could be taken into account satisfactorily. 
Measurements on a long blunt cone (Porter, 1993) then showed that non-uniformity in the 
axial load distribution was no impediment to use of the method. 

Now, for the purposes of stress wave analysis, the relatively slender models which are of 
interest here may be treated as one-dimensional, with cross-sectional area varying in the 
axial direction. Strictly, each axial element of the model produces a particular signal at 
the strain gauges, which therefore give an output which depends on the force distribution 
on the model, as well as the total force. Fortunately, numerical analysis (Simmons et al 
1993), and the experiments noted in the previous paragraph, indicate that for these conical 
models and the times of interest here, the distribution of force is unimportant. When this 
is so, the relation between the output of the strain gauges, y(t), and the force applied to 
the model, u(t), may be written as 

y(t) = u(x) dx , 

Jo 

where g(t) is an impulse response function which relates y(t) and u(t). With g(t) 
determined, either numerically or by calibration, this expression allows numerical 
deconvolution procedures to be applied to y(t) to yield u(t). 

The model employed here was larger and more complicated than the cones of the 
previous studies, so it was necessary to again check the distribution of force was not 
important. This was done numerically. A two-dimensional axisymmetric finite element 
70 x 6 mesh was made to represent the model and fuel tank, and three loading 
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distributions were applied independently to approximately represent the loading on the 
model during a test. One represented the drag of the model, and was applied as a load 
distributed over the surface of the model. One represented the thrust due to injection of 
the fuel, and was applied as a point load at the location of the injection orifices. The 
third one represented the thrust due to mixing and combustion, and was applied as a 
distributed load on the thrust surface. The variation with time of each of these is shown 
in fig. 7(a). Together these three inputs produced a simulated strain gauge response. 
This was then deconvolved, using an impulse response function g(t) which was 
determined for a single point load applied at the nose, to see if the applied time history of 
the force on the model could be recovered. The result is shown as the solid line in fig. 
7(b), where it is compared with the sum of the three forces acting on the model. It can 
be seen that the deconvolved signal is in good agreement with the input signal, signifying 
that the measured force was indeed independent of its distribution. 

Having thus established that no significant error was involved in taking the force as 
applied at the nose, the impulse response function for such a force was determined 
experimentally. This was done by vertically suspending the model and fuel tank, attached 
to the stress bar, by a fine wire attached to the nose. The wire was cut close to the nose 
tip to effect a sudden removal of the tensile load, thereby producing the equivalent of a 
step increase in drag. The resulting strain gauge output shown in fig. 8(a), gave the step 
response for the system, and the impulse response was found by differentiating this with 
respect to time. This approach had the benefit that any effects of stress wave damping, 
though small, would be captured in the impulse response function. With a suitable g(t) 
thus defined in eqn (1), the deconvolution could then proceed numerically, by the time- 
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domain iterative method of Prost & Goutte (1984), to yield the axial force. This process 
yielded a random error of ±10% in the axial force, but no detectable systematic error. 

5. FORCE MEASUREMENTS 

A strain gauge output obtained during a test is shown in fig. 8(b). Fuel injection is timed 
so that a nearly constant fuel mass flow is maintained during the test time, and this has 
the consequence that the thrust due to fuel injection grows before the test time begins. 
This is evident in the lower frequency oscillation in the strain gauge output, before the 
initiation of test flow. As the test flow begins, the frequency increases. 

Using the impulse response function derived from fig. 8(a), the output of the strain gauge 
in fig 8(b) is deconvolved numerical to yield axial force records, as shown in fig. 9(b). 
There are some spurious low amplitude oscillations on these records, which are ascribed 
partly to non-axial vibrations transmitted to the stress wave bar from the surroundings, 
and partly to imperfections in the deconvolution process, but these have only a small 
effect on the accuracy of the records. The test section Pitot pressure, as well as the fuel 
supply pressure for the fuel-on cases, are shown in fig. 9(a), and both are seen to be 
nearly constant during the test time. 

It can be seen in fig. 9(b) that, for the fuel-off case, a drag of approximately 150 N is 
experienced by the model during the test time. When fuel is injected with nitrogen as the 
test gas, the drag is reduced and, before the test flow arrives, a thrust due to the reaction 
of the fuel jets is evident. This thrust increment appears to increase when the flow 
arrives but this is a spurious effect, due to the change in drag associated with the change 
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in test gas. When fuel is injected with air as the test gas, the axial force variation follows 
that of fuel injected into nitrogen, until the fuel ignites in the combustion chambers, when 
the axial force increases in the thrust direction, reaching a plateau at approximately 50 N 
of thrust. 

Measurements of the axial force with and without fuel injection were taken for a range of 
stagnation enthalpies. However, before considering the results of these measurements, a 
theoretical treatment of the effects giving rise to this force is required, and this will be 
done after an analysis of the mode of metering the fuel injection. 

6. FUEL INJECTION 

Initially, it was intended that hydrogen be used as the fuel, as previous experiments 
(Stalker & Morgan, 1984) had suggested that at a combustion chamber length of 60 mm, 
and the pressures of approximately 1 atm, which it was thought would be achieved, 
would be adequate for combustion. However, hydrogen did not bum, due probably to the 
occurrence of expansion waves in the combustor, and a mode of fuel injection which was 
different to the original experiments. Therefore silane was added to the hydrogen to 
promote combustion. 

A system for mixing silane and hydrogen was available, but required excessive time to 
produce a uniform mixture. Therefore the concentration of silane was monitored for each 
test by using the pressure transducer in the fuel supply duct and the pressure in the fuel 
tank, measured before and after a test, together with the following analysis. 


13 


Assuming the flow through the fuel injection orifices is sonic, the fuel mass flow rate at 
any instant may be written as 

m F - {F(i)AJRJ} JM P F = B-]UP f 
where B = F(y) A^RoT, F(y) is a function of the ratio of specific heats, A* is the 

effective area of the orifices, Ro is the universal gas constant, T is the temperature in the 

fuel injection system, which is taken as being constant, M is the molecular weight of the 

fuel, and P F is the pressure measured in the fuel supply duct. From this, the mass Am 

leaving the fuel tank in time t is 

Am = BjM J‘ P F dt . 

Using the change in pressure in the fuel tank, AP r , Am may also be written as 

Am = M.A P r m o T , 

where V is the volume of the fuel tank. Equating the expressions for Am leads to 

M - /;/>,«*/ A?,}’ . 

The fuel supply valve opened in approximately 6 millisec, and closed approximately 
100 millisec later. By monitoring P F during this time, the value of the integral in the 
above equation could be obtained, and AP r could be measured by observing the pressure 
in the fuel tank before and after a test. A previous calibration for B, obtained by 
operating the fuel system with various known amounts of silane in the gas mixture, then 
allowed M to be determined, and hence the silane concentration as well as the mass flow 
through the fuel orifices is also known. This was finally corrected for adiabatic 
expansion in the fuel tank by multiplying the mass flow by 1.08. The equivalence ratio 
was then calculated as the ratio of this mass flow to the mass flow of fuel required to 
consume all the oxygen in the air passing through the duct. 
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7. 


ANALYSIS OF PROPULSIVE PERFORMANCE 


The installed propulsive performance of engines in an aerodynamic vehicle is commonly 
assessed in terms of the thrust generated versus the drag to be overcome. When the 
function of the engine is largely independent of the vehicle aerodynamics, it is possible to 
adopt a simple view, taking thrust as an engine property, and drag as a property of the 
vehicle. However, for scramjets, which are integrated into the flow field about the body, 
such an approach is not possible. In this case it is more suitable to relate the performance 
to fuel injection, estimating the change in axial force on the configuration produced by 
injecting the fuel. 

Therefore 'the analysis proceeds by first considering the drag in the absence of fuel 
injection, before going on to treat the effect of the combustion of the fuel. It should be 
noted that, not withstanding the attempt to use a simple model configuration, the details 
of the flow are so complicated that a rigorous analysis would be extremely difficult. 
Thus, some simplifying approximations are made but, as will be seen below, the analysis 
which results adequately predicts the experimental results. 

8. INVISCID FUEL-OFF DRAG 

Intake Assembly. The flow over the conical forebody and the intakes is shown 
schematically in fig. 10. The conical shock originating at the tip of the cone causes a 
uniform pressure over the region a of the cone surface. The splitters are machined so 
that the inlet ramp surface is always normal to the cone surface, as shown in section AA. 

It is assumed that this leads to shocks which are also normal to the cone surface, and that 
the pressure in the regions b is uniform and is given by the pressure downstream of a 
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two-dimensional oblique shock formed by the leading edge of the inlet ramp. The effect 
of shock-boundary layer interactions is ignored, since they do not affect the overall 
pressure rise, and the dispersion of the pressure rise at the foot of the shock which they 
do cause is likely to raise the force on the surface upstream of the shock by as much as it 
reduces it downstream. As shown in section AA, the outer surface area of the inlet duct 
is greater than the inner area, and this effect is taken into account. 

Calculations indicate that the conical shock from the tip of the cone should pass inside the 
cowl, as shown in figure 10(a). Further calculations showed that if it passed outside of 
the cowl, the leading edge bevel would induce shock detachment. Detachment was not 
evident in- luminosity photographs of the flow, and so the drag could be estimated by 
taking the free stream to be directly incident on the exterior surfaces of the cowl. The 
pressure on the bevel was obtained by treating it as part of an infinite, swept, leading 
edge. 

The maximum temperature in the inlet was less than 2200 K, implying that no 
dissociation took place. Thus the flow model outlined was used, with a ratio of specific 
heats of 1.4. The inviscid drag of the intake assembly, comprising the conical forebody, 
the inlets and the cowl exterior, expressed in Newtons, was calculated to be 

25 x 10r 5 P r 

with the Pitot pressure P p in units of N m' 2 . This was independent of the expected 
variation in the freestream Mach number. The major contributions were the drag on the 
inlet ramps, or splitters, at 50%, the drag on the forecone, at 20%, and the drag on the 
cowl leading edge bevel, at 25%. 
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Combustion Chambers. The inviscid drag of the combustion chambers is zero, since 
their cross-sectional area is constant. However, the pressure and Mach number in the 
combustion chambers is required in order to calculate the level of thrust produced in the 
expansion nozzle. These may be obtained from the relations for adiabatic flow in a 
streamtube (e.g. Hall 1951) which can be written as 

♦iV2(r-i) 


and 


A, M m 

exp (-AS//?) = ~ 

M, 









A . M. 


I 1+1^- Ml 
2 

l + Ml 


\'A 


* 


where A,*, is the cross-sectional area in the freestream of the streamtube captured by the 
intake, M*. and P. are the freestream Mach number and pressure, A f is the total cross- 
sectional area of the combustion chambers, and M f and P f are the Mach number and 
pressure there. Since the conical shock from the cone passes inside the cowl leading 
edges, no "spilling" of the flow from the intake takes place, and A*, can be taken as the 
area defined by the projection of the cowl leading edges onto a plane normal to the flow 
direction. The entropy rise per unit mass of gas associated with the three shocks through 
which the flow is compressed is AS, and R is the gas constant per unit mass. When these 
relations are combined with the dimensions of the model, it is found that, with y = 1.4, 
M f = 4.0 at M„ = 6.1 and increases by 10% for each 10% increase in M a . P f may be 
written as 


P f = (0.56 - 0.059 M_)P p (1) 

which is a form convenient for later use. Estimates indicate that the effect of the 
displacement thickness of the boundary layer on the combustion chamber cross-section is 
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small enough to be ignored. 


Predicted values of P f are compared with measurements in fig. 15. The measurements 
are corrupted by the persistence of waves in the flow, due to intake-combustion chamber 
mismatch, but it can be seen that the predictions are not inconsistent with the mean level. 

Expansion Nozzles. The surfaces of the expansion nozzles are formed by the conical 
afterbody and adjacent splitters, and the pressures experienced on these surfaces 
determine all but a few per cent of the thrust acting on the model. To calculate these 
pressures, it is assumed that the transverse curvature of the conical afterbody may be 
neglected,- and the flow field at the corner between the afterbody and the splitters is that 
produced by double expansion of a streamwise comer, as shown in fig. 11. 

The trace of the wave pattern in two cross flow planes, X and Y, is shown in the figure, 
demonstrating the conical nature of the flow. Taking the plane X, region f corresponds to 
uniform flow which is outside of the zone of influence of the two surfaces. Outside the 
zone of influence of the streamwise comer, the flow passes through one of two centred 
Prandtl-Meyer expansions to produce the two uniform flow regions f, and f 2 . For the part 
of the flow field which is within the zone of influence of the streamwise comer, an 
analytic solution is not known, but numerical solutions (Anderson and Nangia 1977) 
indicate that each of the flow regions f, and f 2 passes through one of two further Prandtl- 
Meyer expansions to produce another region of uniform flow f l2 . The total deflection 
angle experienced in passing from f to f.2 is the same, irrespective of whether the passage 
is accomplished through f, or f 2 . Since the properties of the flow produced by a Prandtl- 
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Meyer expansion are uniquely related to the deflection angle, the properties in f 12 are 
uniquely determined. Thus the flow over each of the two adjacent surfaces can be 
analysed as a uniform flow, produced by the deflection of the surface itself, with the 
deflection of the adjacent surface producing a simple Prandtl-Meyer expansion of that 
uniform flow. 

This model was used to obtain the pressure distribution on the conical afterbody and 
adjacent splitter surfaces, taking account of the expansion wave from the cowl trailing 
edge, and reflection from the plane of symmetry between two splitters, in the same way. 
A simple graphical procedure allowed the pressure distribution to be established with 
sufficient - accuracy. The thrust was then obtained by integration of this pressure 
distribution, taking account of the reduction of area with downstream distance of the 
conical afterbody. The calculations were performed for a perfect gas with a ratio of 
specific heats of 1.4. 

Results are presented in fig. 12 for a range of combustion chamber Mach numbers. For 
comparison, the thrust which would be obtained with one-dimensional flow through a 
nozzle of the same area ratio also is shown, and it is seen that the difference between the 
two increases with Mach number, the thrust of the three-dimensional nozzle falling as the 
one-dimensional one increases. Clearly, the three-dimensional flow makes better use of 
the thrust surfaces at lower Mach numbers than at the higher ones. 

9. VISCOUS FUEL-OFF DRAG 

To assess the viscous drag of a surface, it is first necessary to determine whether the 
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boundary layer is laminar or turbulent. Previous experiments (He and Morgan 1994) 
have shown that transition occurs on a flat plate in the shock tunnel at a Reynolds’ 
number, based on distance from the leading edge, of approximately 2 x 10 s . Applying 
this to the model used in the experiments, the boundary layer flow on the inlet splitter 
face, the internal inlet, surfaces of the cowl, and the exterior surface of the cowl, along 
its full length, is expected to be laminar. The Reynolds’ number on the conical forebody 
may reach 4 x 10 6 but, allowing for the higher transition Reynolds’ numbers on a cone, 
this too is expected to exhibit predominantly laminar flow. 

Thus for these surfaces, the local skin friction coefficient is given by (e.g. Hayes and 
Probstein 1959) 

9 = 

- 0.664 n(p' u ' / P. ! 

where p e , u, and p c are the density, velocity and coefficient of viscosity in the external 
stream adjacent to the boundary layer, r w is the skin friction per unit area, is the 
Reynold’s number based on the external stream properties and the distance from the 
leading edge, n = 1 and V3 for a flat plate and a cone respectively, and p* and p* are the 
density and viscosity at a reference temperature T\ given by 

T'/T e = 0.5(1 + TJT e ) + 0.0384 

where T e and T w are temperature in the external stream and at the wall respectively, and 
M e is the external stream Mach number. Due to the reduction in Reynolds’s number as 
the stagnation enthalpy increases, the skin friction increases with stagnation enthalpy, H,, 
and the total drag, in Newtons, for the surfaces noted in the paragraph above is 

(3.09 + 0.55 HJ x ia 5 P r 
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where H, is in MJ/kg. Of this the forecone, the inlet splitter faces, the interior inlet 
surface of the cowl, and the exterior surface of the cowl til contributed approximately 
equal portions. 


The Reynolds’ number of the flow in the combustion chambers is such that, taking 
account of the shock-boundary interactions in the inlets, the boundary layer flow is 
expected to be turbulent. A review of turbulent skin friction by Hopkins and Inouye 
(1971) indicates that the theory of Spalding and Chi is most suited to the low values of 
T„/T e experienced here. This theory obtains the skin friction in compressible flow from 
incompressible values by using the transformations 

}f = f . c /> 
and R a = F s R a , 

where the bar denotes the incompressible values. F c and F, are given here by 

F c = 0.18 M 2 ' (sin' 1 a + sin' 1 ?)' 2 , 
andF x = (sin' 1 a + sin l pf {O.IS M? (TJT e f ™ (TJTJ) 0 ™}' 1 
where a = (2/1 2 - B) (44 2 + B 2 )~° s , 
and p = B{4A 2 + B 2 )^" 5 , 
with A 2 = 0.18 U) TJT W , 

and B = (l + 0.18 Af 2 - TJT e )TJT w . 


T w is the stagnation temperature in the external stream. The incompressible values of 
skin friction were taken from the relation 

c, = 0.060 R* 2 . 

The effective origin of the turbulent boundary layer was obtained by estimating the 
laminar boundary layer thickness at the entrance to the combustion chambers after 
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compression to the combustion chamber pressure, using this to indicate the effective 
origin of a laminar boundary layer, *nd assuming that the turbulent boundary layer 
originated at half that distance upstream. Using this theory, the skin friction drag in the 
combustion chambers, expressed in Newtons, was calculated to be 

(35.1 + 0.82 HJx 10 s P f 

where P f is in units of N m' 2 . There is negligible variation over the range of Mach 
numbers experienced without fuel injection. 

The boundary layer flow in the expansion nozzles was also taken to be turbulent, as the 
favourable pressure associated with the expansion comer was estimated as insufficient to 
cause relaminarisation (Narasimha and Viswanath 1975). Noting that the skin friction 
coefficient varies slowly with Reynold’s number it is assumed to be unchanged across the 
expansions and therefore, by also assuming that the velocity is unchanged, the skin 
friction at any point on the surface of the expansion nozzles can be obtained by knowing 
the local density in the inviscid flow. This can be derived from the pressure distribution 
on those surfaces. Integration over the surfaces then yields the drag due to skin friction 
on the expansion nozzles,in Newtons, as 

4.8 x ia 5 P f . 

The skin friction on the part of the model support sting exposed to the flow was estimated 
to be negligible. 

10. RESULTANT FUEL-OFF DRAG 

Fig. 13 displays the results of the above calculations, with eq(l) used to relate P f to P p . 
Fig. 13(a) shows the inviscid drag, as the drag due to the intake and the cowl, minus the 
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thrust due to the expansion nozzle. Ideally, it should be possible to design a model such 
that one largely canals the other, apart from a relatively small effect due to entropy rise 
in the intake, and to the reduction in expansion nozzle area ratio occasioned by the need 
to support the model on the sting. The fact that the residual drag is not small is largely 
due to two effects. The first is the drag due to the cowl leading edge bevel, which was 
incorporated solely for ease of manufacture of the model, and the second is the limited 
thrust performance of the expansion nozzle. This is due partly because it is operating at 
combustor flow Mach numbers about 4, which, as shown by fig. 12, is considerably in 
excess of the optimum value, and this entails a considerable reduction in thrust. It is also 
due to the fact that even at lower Mach numbers, its thrust performance could be 
improved, -as shown by the comparison with a one-dimensional nozzle in fig. 12. 

Fig. 13(b) shows the viscous drag. It can be seen that the drag divides approximately 
evenly between the combustion chambers, where the boundary layer is turbulent, and the 
intake plus the cowl exterior surface, where the boundary layer is laminar. The 
expansion nozzle makes a small contribution only. The total viscous drag is comparable 
with the inviscid drag which, as already mentioned, can be substantially reduced by good 
design to leave viscous drag as the dominant contribution. This illustrates the point that 
successful operation of scramjet vehicles will be heavily dependent on control of these 
viscous losses. 

The inviscid drag and the viscous drag are combined to yield the resultant drag in fig. 
13(c), where it is compared with measurements. The possible error in the measurements 
is about equally due to uncertainty in the measured drag and uncertainty in the measured 
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Pitot pressure. Within the limits of error, it is clear that theory and experiment agree, 
except for a tendency for the drag to increase somewhat more rapidly with stagnation 
enthalpy than is indicated by the theory. This could be due to the theoretical model used 
for turbulent skin friction. 

11. INVISCID FUEL-ON THRUST 

Combustion chamber pressure rise. The propulsive effect of a scramjet derives from 
the pressure rise experienced in the combustion chambers due to mixing and burning of 
the fuel there. If either, or both, of these processes are not completed in the combustion 
chamber the pressure rise is reduced, and therefore so is the thrust. 

Shock tube data (Jackimowski & McLain 1983) indicates that, regardless of the presence 
of silane, the fuel should not bum in the combustion chambers. However shock tunnel 
experiments (Morris 1988) exhibit much shorter reaction lengths, due possibly to high 
boundary layer temperatures and "frozen" oxygen radicals in the freestream. Under the 
conditions of the present tests, reaction lengths of the order of 10 mm are likely, and as 
this is much less than the length of the combustion chambers it is assumed that complete 
combustion of all the fuel that is mixed takes place. 

An indication that essentially complete mixing occurred was obtained by plotting the 
thrust obtained against the fuel equivalence ratio, as in fig. 14(a). It can be seen that the 
thrust varied approximately as the fuel mass flow, suggesting that the proportion of fuel 
burnt was approximately independent of the mass flow. Since the proportion of a fuel jet 
which is mixed in the length of combustion chamber will only be invariant with the jet 
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mass flow if all of the jet is mixed, it is reasonable to assume that complete mixing 
occurred. 

Thus, an analysis was conducted of pre-mixed, equilibrium combustion in a constant area 
duct of the silane-hydrogen fuel. It was assumed that the silane burned according to the 
reaction 

SiH A + 20 1 - Si0 2 + 2H 2 0 , 
and the hydrogen according to the reaction 

2 H 2 + 0 2 - 2 H 2 0 . 

Dissociation of H 2 0 in the combustion products was taken into account (e.g. Beniot 1968) 
and, since post combustion temperatures were above the condensation temperature of 
Si0 2 , it was assumed to remain in the vapour state. The analysis was conducted in a 
conventional manner, using the energy equation, with the heats of formation included in 
the enthalpy, together with the momentum equation and the equation of state, applied 
across a shock-like discontinuity. The change in molecular weight due to addition of the 
fuel, and due to its chemical reaction with air, were taken into account. For each set of 
precombustion conditions, a solution was obtained iteratively, yielding the results shown 
in fig. 14(b). They show that the combustion heat release has a reducing effect as the 
stagnation enthalpy, and therefore the pre-combustion temperature, is increased. The 
combustion pressure rise is reduced, and the post-combustion Mach number increases, to 
come closer to pre-combustion values. A curve for the temperature ratio across the 
combustion zone is also presented, showing that the pressure rise is almost entirely due to 
the combustion temperature rise. 
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The effect of dissociation of H 2 0 in the combustion products on the temperature rise 
increased with stagnation enthalpy, reducing T c by 6% at 3 MJ kg ', and 25% at 
8.5 MJ kg' 1 . Thus, at the higher stagnation enthalpy, dissociation reduced P^Pf from 1.5 
to 1.2, which is not a large effect when compared with that of pre-combustion 
temperature. Hence, for the conditions of this paper, dissociation is not considered to be 
important. 

The effect of combustion was measured by making a series of pressure tappings in one of 
the combustion chambers, as shown in fig. 15. Although the non-uniformity of the flow 
makes an exact comparison with the theory impossible, it can be seen that if the average 
of the two* downstream data points is taken, the pressure level is predicted reasonably 
well. The upstream data points may be taken to indicate that mixing and/or combustion 
is not complete at those stations. As shown, Pitot pressure was also measured at a point 
which was just upstream of the expansion wave from the comer of the centrebody. 
Taking the mean of the static pressure at the two downstream stations when combustion 
occurred, and the ratio of specific heats for the combustion products, and combining these 
with the Pitot pressure, it was possible to obtain the Mach number of the flow after 
combustion. This varied from 1.5 to 2.2 as the stagnation enthalpy increased from 
3.2 MJ kg' 1 to 5.8 MJ kg' 1 , signifying that supersonic combustion had indeed taken place. 

Expansion Nozzles. The thrust due to the expansion nozzles is obtained as before with 
P c replacing P f as the pressure before expansion. The frozen ratio of specific heats for 
the combustion products differed from 1.4 by only 3%, so fig. 12 could be used again to 
predict the thrust. However, it should be noted that combustion reduces the Mach 
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number at the entrance to the nozzles, an effect which was taken into account by using 
the theoretically calculated Mach numbers presented in fig. 14. 

Fuel Injection Thrust. There is a substantial amount of thrust associated with fuel 
injection, and to take account of this it was noted that the fuel supply duct pressure, P F , 
was an order of magnitude greater than the pre-combustion pressures in the combustion 
chambers. Thus it was assumed that after injection, the fuel expanded to very low 
pressures. The velocity of the fuel could thus be obtained as ^2/( y-1) a F , where a F is 
the speed of sound in the fuel before injection, and by combining this with the fuel mass 
flow, and allowing for the injection angle, the thrust could be determined. A discharge 
coefficient 'of 0.9 was used for the fuel orifices, and for fuel mass flows corresponding to 
an equivalence ratio of 0.77, the thrust in Newtons was 

(60 - 4.2 HJ x 10 s P r . 

where eq(l) can again be used to express P f in terms of P p . While this was not 
inconsistent with measurements of thrust before the test flow arrived, as seen in fig. 9, 
the accuracy of those measurements prevented a more definitive comparison. 

The results of these calculations are presented in fig. 16(a), for an equivalence ratio 
corresponding to the mean value used in the experiments. The inviscid drag of the intake 
and cowl remains as before, but the pressure rise due to combustion allows the expansion 
nozzle to produce much greater thrust. When this is added to the fuel injection thrust it is 
sufficient to outweigh the inviscid drag over most of the range of enthalpies studied. 
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12. VISCOUS FUEL-ON DRAG 


The injection of fuel does not, of course, influence the skin friction drag on the intake 
and cowl, but it does affect the skin friction in the combustion chambers and the 
expansion nozzles, through the effect of combustion on the temperature of the flow 
external to the boundary layer. The theory of Spalding and Chi is used once again to 
yield 

(70 - 3.0 HJx ia s P f 

for the skin friction drag, in Newtons, on the combustion chambers under the calculated 
flow conditions. The skin friction drag on the expansion nozzles is also estimated 
according to the same approximation as in the fuel-off case, to yield ^in Newtons, 

6.3 x ia s P f . 

These results are presented in fig. 16(b) in terms of P p . Comparison with fig. 13(b) 
reveals that the major change in skin friction drag associated with fuel injection is the 
increase in drag of the combustion chambers. This amounts to a 90% increase at a 
stagnation enthalpy of 3 MJ kg* 1 , where the combustion has the largest proportional effect 
on the temperature. However, at 8 MJ kg' 1 , where the proportional increase in 
temperature is greatly reduced, the increase in combustion chamber skin friction is only 
14%. 


The skin friction in the expansion nozzles is also increased, but this is such a small 
contribution to the overall skin friction drag that the magnitude of the increase can barely 
be noticed on the scale of fig. 16(b). 
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13. RESULTANT FUEL-ON AXIAL FORCE 


The curves in fig. 16 for total inviscid axial foice, and total viscous drag, are combined 
into one curve for the resultant axial force in fig. 17. In order to finally emphasise the 
propulsive performance, the curve is presented with thrust as the positive ordinate. 
Measurements of thrust are also presented, for fuel equivalence ratios near the value used 
in the theoretical calculations, and are seen to follow the theoretical predictions 
reasonably closely, with a tendency to fall below the theoretical curve at the higher values 
of stagnation enthalpy. This follows the same trend as for the fuel-off case, and to 
facilitate the comparison, the resultant drag curve and the measurements of fig. 13(c) are 
presented again in fig. 17. It seems likely that the cause of the discrepancy is the same in 
both areas. * 

It will be noted that the maximum thrust is obtained at the lowest stagnation enthalpy. 
Attempts to increase the maximum value achieved, either by lowering the stagnation 
enthalpy or by increasing the fuel equivalence ratio, produced unsteadiness in the force 
measurement and then "choking" of the flow through the model. In view of the low 
value ( ~ 1.5) measured for the post combustion Mach number at a stagnation enthalpy of 
3.2 MJ kg' 1 , it seems reasonable to conclude that this was due to thermal choking of the 
flow in the combustion chambers. 

The propulsive performance of the model used in these experiments was not outstanding. 
For example, if the maximum value of thrust obtained is taken as 10 x Id 5 P p Newtons, 
the corresponding specific impulse is 175 sec. If the silane in the fuel were replaced by 
an amount of hydrogen sufficient to consume the same amount of oxygen, and the same 
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thrust were obtained, the fuel specific impulse would be 370 sec. However, it is 
worthwhile reiterating that tie model was designed for ease of manufacture and analysis, 
and not for optimum performance. 

Inspection of figs. 13 and 16 suggest ways in which the performance might be improved. 
Fig. 13(a) indicates that the inviscid drag of the configuration could be reduced; for 
example, by removing the external bevel on the cowl leading edges. The thrust delivered 
by the expansion nozzles could be improved by better design, although it must be 
remembered that combustion lowers the Mach number of the flow entering the expansion 
nozzles, and as shown in fig. 12 this improves the efficiency of the nozzles. Also 
expansion* nozzles of larger area ratio would deliver more thrust, though provision of a 
larger expansion area ratio would impact other aspects of the overall design. 

The viscous drag, as shown in figs. 13(b) and 16(b), represents a substantial portion of 
the overall drag, and one which offers less opportunity for reduction. However, there is 
a clear benefit in this respect in reducing the length of the combustion chambers, provided 
this can be done without reducing the benefits of combustion. 

The relatively large size of the viscous drag points to one of the major difficulties 
associated with air breathing propulsion at very high speeds. It is likely to be an 
intractable form of drag, difficult to design against, and its importance here confirms the 
widely held view that it may be a major impediment to effective propulsion at such 
speeds. 


Figs. 16 and 17 display the rapid decrease in net thrust which is occasioned by increasing 
stagnatio.i enthalpy. This is, of course, due to the increasing combustor temperature 
before combustion, which reduces the combustion pressure rise. For application to an 
acceleration vehicle, a scramjet must operate near its optimum over a range of stagnation 
enthalpies, and therefore a means of maintaining the temperature before combustion 
reasonably constant would need to be found. 

14. CONCLUSION 

Thrust and drag measurements have been made in a shock tunnel on a supersonic 
combustion ramjet model. The effect of stagnation enthalpy on the propulsive 

performance of the model was measured for a range of stagnation enthalpies 
corresponding to flight velocities up to approximately 4 km/sec. A stress wave force 
balance was used, which involved monitoring the stress waves produced by the forces on 
this model as they passed into the model sting support. This, in combination with a 
pulsed fuel injection system, made the measurements possible. 

The results were compared with theoretical predictions, and were found to be in generally 
good agreement, with a tendency to slightly underestimate the drag at high enthalpies. 
The analysis leading to the predictions involved the use of established theories which, 
together with some reasonable simplifying assumptions, provided estimates of the relative 
importance of the factors affecting the thrust and drag. Except for the combustion 
process itself, the analysis assumed a perfect gas with a ratio of specific heats of 1.4. 

The experimental validation of this relatively simple theoretical analysis may be expected 
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to have a number of consequences. As well as indicating the significance of factors such 
as expansion nozzle performance and skin friction drag, it allows the benefits of concepts 
for improving the propulsive performance to be quantified, and thereby assists in 
improving the performance. Also, as with any flight vehicle, a simple analysis based on 
results of model experiments in a wind tunnel is useful for preliminary estimates of 
performance. Thus, with proper attention to factors such as the location of transition, 
estimates may be made of the performance of a scramjet powered vehicle in flight. 

Using the theory, it was possible to identify areas where improvement may be expected to 
yield a gain in the overall thrust. For acceleration vehicles, it is clearly of primary 
importance that a means should be found of maintaining the precombustion temperature 
constant as the flight velocity changes. For cruise propulsion, this is less important. In 
general the inviscid drag of the configuration could be reduced by attention to design 
details and the performance of the expansion nozzles could be enhanced by matching them 
to the Mach number at which they are expected to operate. The theory also revealed the 
importance of viscous drag, and drew attention to the value of reducing viscous drag by 
reducing combustion chamber length. 

However, perhaps the most important result is that the investigation has established a 
method for experimental testing of scramjet powered models at high velocities, together 
with the theoretical understanding which attends the approximate analysis. These two 
factors together may be expected to assist in the development of scramjet research. 
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UST OF FIGURES 


Fig 1 . Free Piston Shock Tunnel T4 

Fig 2. Shock tunnel test time 

(a) Typical nozzle reservoir pressure record,!!, = 7.7 MJ kg' 1 

(b) Mass spectrometer measurement of driver gas contamination. 

te = time after start of flow in test section until 10% driver gas 

contamination 

H, = stagnation enthalpy. 

Fig 3. Test Section Flow Conditions 

(a) Measured Pitot pressure 

(b) Calculated Conditions. 

P p = Pitot pressure, M w = Mach number, T*, = static temperature, 
u* = velocity. 

Fig 4. Test Section Flow Calibration and Location of Model. Stagnation 
enthalpy = 8 MJ kg 1 
• Pitot survey measurements 

d = downstream distance from nozzle exit, r = distance from centreline. 

Fig 5. The Scramjet Model 

Fig 6. Model, Fuel Supply System and Model Mount. 

Fig 7. Simulation of Distributed Load 

(a) Input loads. fuel injection Drag Thrust. 

(b) Output. deconvolved load true input. 

Fig 8. Strain Gauge Outputs 

(a) Response to step load 

(b) Response during test. 

Fig 9. Axial Force Measurements. Stagnation enthalpy = 3.2 MJ kg' 1 . 

Fuel-on Air equivalence ratio = 0.83 

(a) Test flow and fuel monitoring 

(b) Force on model. 

Fig 10. Intake Flow (i) Forecone and inlet cowl. 

(ii) Shock pattern between splitters. 

Fig 11. Expansion Nozzle Flow 

Fig 12. Thrust Nozzle Performance. Area ratio = 5.84, y - 1.4 

P f = combustion chamber pressure, M f = combustion chamber Mach number. 
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Fig 13. 


Fig 14. 


Fig 15. 


Fig 16. 


Fig 17. 


Drag Without Fuel Injection. theory, a experiment. 

X 

H, = stagnation enthalpy 

(a) Inviscid Drag 

(b) Viscous Drag 

(c) Resultant Drag. 

Combustion Features 

(a) Experimental effect of fuel mass flow. 

(b) Theoretical effect of combustion. 

AF = thrust increment due to fuel injection, <f> ~ equivalence ratio, P c = post 
combustion pressure, T c = post combustion temperature, M c = post 
combustion Mach number, H, = stagnation enthalpy. 

Combustion Chamber Pressures 

Fuel-off: O Experiment, theory. 

Fuel-on: a Experiment, theory (<t> = 0.77) 

1 

H, = stagnation enthalpy. 

Axial Force with Fuel Injection and Combustion 
Equivalence ratio = 0.77, F A = axial force for thrust), 

H, = stagnation enthalpy 

(a) Inviscid axial force 

(b) Viscous Drag 

Resultant Axial Force 

Fuel-on: Equivalence ratio <t> - 0.77 ± 0.10 

V Experiment, theory. 

X 

Fuel-off: a Experiment, — theory. 

X 

F a = axial force, P p = Pitot pressure, H, = stagnation theory. 
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